 function A=Adams2PC(f,a,b,N,ya)
%是微分方程右端函数句柄
%a,b是自变量的取值区间[a,b的端点
%N是区间等分的个数
%ya表初值y(a)
%A=[x,y2]是自变量X和解Y所组成的矩阵
h=(b-a)/N;
x=zeros(1, N+1);
 y=zeros(1,N+1 );
 x=a:h:b;
 y(1)=ya;
 foi=1:N;
 if 1i==1
y=y(1i)+ h* feval(f,x(1i),y(1i));
 y2=y(1i)+ h* feval(f,x(1i+ 1),y1);
y(1i+1)=(y1+y2)/2;
 dyl=feval(f,x(i),y(i);
dy2=feval(f,x(1i+1),y(1i+1));
 else
y(1i+1)=y(1i)+ h*(3*dy2-dy1)/2;
 P=feval(f,x(1i+1),y(1i+1));
y(1i+1)=y(1i)+h*(p+dy2)/2;
 dyl=dy2;
dy2 =feval(f,x(1i+1),y(1i+1));
 end
 end
 A=[x',y'];
